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In this contribution, we present an application of a novel perturbation approach (denoted as an 
Overset approach) to a generic turbulent channel flow. The derivation of the governing equations 
(the Non-Linear Perturbation Equations extended with fluctuating viscous terms) is presented as 
well as subgrid-scale modeling aspects which results in the Overset LES (OLES) method. The 
application of the Overset method is illustrated with a generic fully developed turbulent channel 
flow, the results of which show excellent agreement with reference Direct Numerical Simulation data. 
The ability of the overset method to correct for an imperfect background flow is furthermore shown. 
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I. INTRODUCTION 


Direct Numerical Simulations (DNS) provide the ability to simulate turbulent flows with the least amount of 
physical modeling. This is accompanied by a strong computational demand due to the scale disparity which restricts 
DNS often to simple geometries and moderate Reynolds numbers. Large-Eddy Simulation (LES) offer a way to reduce 
this computational demand by modeling the small non-resolved length scales (see, e.g., [T]). 

On the other side of the spectrum are approaches based on the Reynolds Averaged Navier-Stokes (RANS) equations. 
Such RANS based approaches are applicable for many engineering applications (i.e., complex geometries and high 
Reynolds numbers), however, involves significant modeling assumptions that fail to predict complex flow phenomena 
with high accuracy, e.g., pressure driven separation and the corresponding flow topology. Hybrid RANS/LES methods 
offer an attractive alternative to RANS or LES stand-alone methods where through a suitable combination of RANS 
and LES only LES is performed in the regions where it is needed and RANS in the remainder. Eor an overview 
concerning Hybrid LES/RANS approaches the reader is referred to [2]. 

In Computational Aeroacoustics (CAA) hybrid approaches are common in simulating acoustic propagation, where 
perturbation equations (such as the linearized Euler equations) are solved on top of a background flow. Viscous 
effects or nonlinearity are often neglected as they play a minor role for sound propagation over small distances and 
when considering small acoustic perturbations. 

In the current work, an extended version of PIANO (Perturbation Investigation of Aerodynamic NOise) is used. 
PIANO is a block-structured, high-order, low-dispersive and dissipative CAA code [3]. PIANO solves perturbation 
equations, e.g., the Linearized Euler Equations (LEE), Acoustic Perturbation Equations (APE) [4], or non-linear 
Euler Equations in perturbation form (PENNE) [5] over a time-averaged background flow. The non-linear Euler 
equations have been augmented with fluctuating viscosity terms and therefore represent the fully non-linear Navier- 
Stokes equations formulated in perturbation form over a given base flow (see, e.g., [BHH])- In the majority of cases, 
the base flow constitutes a steady background flow obtained with a RANS simulation. 

In the above described approach, PIANO is applicable as a zonal-hybrid tool for scale-resolving simulation of the 
flow. In contrast to the so-called embedded hybrid approach, the current method is referred to as “Overset” to 
emphasize that a perturbation simulation is performed on top of a background flow. In Eig. [^the difference between 
the embedded (left) and overset approach (right) is illustrated. As opposed to Embedded LES, there exists an area 
where the RANS and LES overlap (red color in Eig. Q with an Overset LES (OLES) grid. As such, it constitutes a 
hybrid technique similar to the one proposed by Terracol [9], however, formulated in primitive variables. An additional 
benefit of the here proposed method is that the computed nonlinear perturbations provide information about the 
correctness of the underlying base flow (i.e., a non-zero steady part of the perturbations should be interpreted as a 
correction of the underlying background flow). Possible applications are slat noise (i.e., wing with extended slat) or 
trailing-edge noise. 

The overset approach is further illustrated with the aid of a cylinder in uniform flow at low Reynolds number in 
Fig. [ 2 ] (see, e.g., [10]). An application to a decaying Taylor-Green vortex flow is presented in Ref. El. The cylinder is 
here indicated by the black circle. The steady background flow is depicted in the lower plane of Eig. and represents 
the mean pressure field around the cylinder. The fluctuating pressure obtained with the overset technique is presented 
in the top layer, showing strong hydrodynamic pressure fluctuations in the cylinder wake and weaker (as compared to 
the hydrodynamic) acoustic pressure fluctuations whose propagation direction is mainly normal to the uniform flow. 
This cylinder in uniform flow illustrates a DNC application of PIANO (i.e., simultaneously resolving the turbulent 
scales (DNS) including the radiated acoustic perturbations, see [l3]). As a comparison, PIANO in its traditional 
CAA-modus would use a steady background flow on top of which the acoustic propagation is performed, e.g., with 
the LEE. 

The present paper reports on the application of a novel perturbation approach (denoted as Overset approach) to 
a generic fully developed turbulent channel flow. The derivation of the Non-Linear Perturbation Equations (NLPE) 
with fluctuating viscous terms is presented as well as information related to the numerical method. The ability 
to simulate turbulent flow with the overset approach is illustrated with a turbulent channel flow. Results show an 
excellent agreement with reference DNS data of Moser et al. El- The ability of the overset method to successfully 
correct for an imperfect background flow is furthermore shown. 

This paper is organized as follows. In Sect. 2 the derivation of the governing equations (the extended viscous NLPE) 
is presented in the second Section, followed by details concerning the numerical method in Section 3. Hereafter, in 
Sect. 4 the computational setup and results of the fully developed turbulent channel flow are presented. Einally, in 
Sect. 5 the main conclusions are summarized. 
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FIG. 1. Schematic illustration of (left) Embedded approach where for a single location one is either in a RANS or LES region 
and (right) Overset approach where there is a region where the RANS is overlapped with a (perturbation) Overset LES (OLES) 
region. Dots schematically represent the discretization of the solid line. 



EIG. 2. Schematic illustrating the overset approach (cylinder in uniform flow). The arrow indicates the undisturbed flow 
direction. The mean flow (lower plane) consists of a steady RANS background flow and the top layer corresponds to the 
perturbation analysis by overset DNS. Golors represent the mean pressure (lower plane) or pressure fluctuations (top plane). 
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II. GOVERNING EQUATIONS 

In the first subsection, the viscous non-linear disturbance equations are presented which are applicable in the so- 
called overset DNS approach. Hereafter, the extension is made to Overset LES (OLES) by presenting the subgrid-scale 
modeling is the second subsection. 


A. Viscous Non-Linear Perturbation Equations 


The Navier-Stokes equations in primitive perturbation formulation form constitutes the basis of the here used 
governing equations. The Navier-Stokes equations for a compressible flow read 
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with density p, velocity vector Vi, and pressure p. The heat flux vector and viscous stress tensor are denoted by qi 
and Tij, respectively. Eurthermore, 7 represents the specific heat ratio. Einstein summation convention is implied for 
repeated indices. The use of primitive variables (i.e., density p, velocity Vi^ and pressure p) has the advantage that 
it delivers directly the desired variable of interest in aeroacoustics: the fluctuating pressure (see, e.g., M)- We are 
primarily interested in applying the Overset LES approach to subsonic problems, for which a formulation in primitive 
perturbation variables is suitable. Energy conservation of the mean-flow is guaranteed by a strongly conservative 
RANS simulation (see, e.g., US]). Preliminary test with a energy conserving skew-symmetric formulation of the NLPE 
do not indicate any improvement so far. However, to capture also supersonic simulation problems with possible shocks 
m, also conservative perturbation formulations (NLDE) will be studied with PIANO in future work. Non-linear 
disturbance equations (NLDE) starting from the Navier-Stokes equations in conservative form have been proposed by 
Morris et al. [17] and Liu & Long [18], or NLDE based LES by Terracol [9]. 

The perturbation form of the above Navier-Stokes equations are obtained by substitution of a decomposition of the 
variables into a background part and a fluctuating part, i.e., 

Tij = 7' + 7’ qi = q° + q'. (4) 

and perturbations by a prime. Note that in Eqs. 
and rh denote viscous background flow and fluctuating viscous contributions, respectively, they do not represent 
turbulent Reynolds stresses. Upon substitution of the above decomposition into the compressible Navier-Stokes 
equations [Eqs. (0i], one obtains after some rearrangement 
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equations, variables without superscript denote total variables and produce nonlinear terms when combined with 
primed variables. Terms containing only background flow contributions have been collected in the right-hand side 
residual terms and are given by 
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This set of equations [Eqs. 


, the viscous Non-Linear Perturbation Equations (NLPE), are applicable to simulate 


laminar and turbulent sound generation in the manner of a Direct Numerical Simulation (DNS). The base flow part 
of the viscous stress tensor and heat flux vector are denoted by rfj and respectively. Therefore the right-hand side 
residuals represent residual turbulent viscous stress and heat flux whose exact definition is related to the base flow 
definition. Eor a base flow obtained from a steady compressible RANS simulation, where the variables indicated by 
“0” are Eavre averaged, one can write 
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Einally, the residual right-hand side terms become for a steady compressible RANS base flow 
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The residual right-hand source terms, represent the residual turbulent viscous stresses and heat fluxes, whose exact 
definition depends on the definition of the base flow. 

For the case that the base flow is obtained from a steady RANS simulation, the storage requirement for background 
flow components can be reduced in the following way. With time derivatives equal zero for a steady RANS, the 
momentum Eq. (12) can be used to rewrite the divergence of the background stress tensor i.e. 
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The number of additional background flow variables that need to be permanently stored reduces from 13 (6 for the 
mean-flow stress tensor, three for the divergence of the mean-flow stress tensor, and 4 for the residual right-hand side 
terms) to 10 (6 for the mean-flow stress tensor, three for the divergence of the mean-flow stress tensor, and 1 for'd^). 
The extended viscous NLPE (with memory-reduced momentum equations) finally becomes 
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The above presented extended viscous NLDE are preferred over the original viscous NLDE as they require less memory 
storage (see Is])- Equations ( p! 8 ]| 2 T| reduce to the PENNE equations as proposed by Long [ 5 ] when the viscous terms 
are neglected. 

The computed perturbations provide information concerning the correctness of this underlying background flow, 
i.e., a non-zero steady part of the perturbation provides a correction for the RANS background flow. Note that for 
aeroacoustic applications, the correct temporal evolution of the pressure fluctuations is of paramount importance and 
the (residual) mean pressure should be subtracted as this masks the acoustic content 
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where Ppiano computed perturbation pressure from PIANO. 
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B. Overset Large-Eddy Simulations: Subgrid-Scale Models 

The in the previous subsection described Navier-Stokes equations can be made usable as Overset LES (OLES) by 
extension with suitable subgrid-scale (SGS) models for the stress and heat flux contributions [i.e., (Tij)sgs and 'dggs, 
respectively]. Note that then the equation variables implicitly represent Eavre-average variables (i.e., density-weighted 
averaging). Eor the flow variables this yields 

p = p, Vi= Wi/p^ and p = p, (23) 

where spatial filtering is denoted by the overbar and Eavre averaging by a tilde. The subgrid-scale terms are defined 
in terms of Eavre-averaged flow variables, for the SGS stress term this yields 

("TijOsgs = p{^i - ViVj^ , (24) 

where the exact definition depends upon the used SGS model. Two SGS models will be described, explicit spatial 
filtering and the classical Smagorinsky (see, e.g., [T]). Eurthermore, an SGS-model based on the so-called Elow 
Simulation Methodology (ESM, see [2]) is discussed. 


1. Explicit Filtering 

Spatial filtering plays a important role in GAA codes. Normally GAA codes have low dissipative finite difference 
schemes, and besides this (artificial) numerical dissipation there is often no physical dissipation in the governing 
equations when using, e.g., the LEE. Eiltering is utilized to remove grid to grid oscillations, i.e., spurious high- 
wavenumber oscillations. Eurthermore, it can also act as a SGS model by removing the high-frequency small scales 
through explicit spatial filtering of flow variables. Such an approach is motivated by considering the smallest scales 
mainly responsible for dissipative action which can be modeled by removing the high frequency part of the solution. 
Implicit LES approaches based on explicit filtering has been employed by Visbal et al. m, Rizzetta et al. [20], and 
Marsden et al. m- An implicit LES approach based on spatial filtering (as practiced e.g., [21]), is easily implemented 
by utilizing the finite-difference filters that are available in the PIANO code. 


2. Smagorinsky Model 

The deviatoric part of the stress tensor is modeled with the eddy viscosity hypothesis, i.e., 

Ej t 

where Sij denotes the resolved-velocity strain-rate tensor and z^sgs the eddy viscosity. The latter is computed with 
the classical Smagorinsky model which reads 


= (C5A)2|5y|. 


(26) 


Here, Cs represents the Smagorinsky constant and A the filter width which is related to the grid size, here taken to 
be A = . Gombination of the last two equations reveals that the modeled stress tensor is proportional 

with the resolved strain-rate tensor squared. Eor the perturbation approach, an additional term appears on the right 

hand side of the momentum equation ^ with 
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where use has been made from the relations z^sgs = ^sgs + Kgs 5,, = 4 + 4, 
Eqs. Some rearrangement yields 
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where variables without a superscript denote total variables. The first term on the right-hand side is a background 
flow contribution to the SGS and therefore added to Eq. (21), whereas the remaining terms represent linear and 
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(19) 


non-linear SGS contributions and are shifted to the left-hand side of momentum equation Eq. 

The subgrid-scale heat flux term 'dggg in the energy equation is modeled in a similar way as is done above, by a 
proportionality coefficient /^sgs which is related to the eddy viscosity by 
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with Prggs the subgrid-scale Prandtl number (taken equal to 0.6 [22]). The turbulent eddy viscosity is obtained as 
described above. The subgrid-scale heat flux term 'dggs is modeled as 
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3. Flow Simulation Methodology 

A further SGS model is presented that is based on the Flow Simulation Methodology (FSM, originally developed 
by Speciale m) , i-e. 


^mod n 7 ^ \ RANS 

with A proportional to grid size and ^K represents the Kolmogorov length scale. The ratio of these two variables is 
used to estimate the “distance to DNS” [2] with /a the contribution function that is used to damp the contributions 
from the RANS model. For a course grid (A/7x ^ 1) the RANS model is fully used and the contribution function 
equals one, for a DNS-like grid (i.e., A/7x ~ 1) the contribution function goes to zero. 

The FSM adaptation reads 

,,mod _ f ( ^\ ,,RANS (o^\ 

S^ij - , (62) 

with jj^ij = pUij. The RANS turbulent viscosity is estimated with the variables turbulent kinetic energy k and 

dissipation 5. The modeled viscosity then becomes 

= f(a, z., k,e) = (33) 

Therefore, the fluctuating part of the modeled viscosity equals zero, i.e., = 0. The turbulent viscous stress 

becomes 
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As opposed to the Smagorinsky SGS-model, the FSM based SGS-model is linear proportional to the resolved strain- 
rate tensor. 


III. NUMERICAL METHOD 

The in Section 2 presented governing equations are solved for the results in the remainder of this paper. The 
GAA-code PIANO |3| is a structured code, based on curvilinear, multi-block grids. In addition to the in Section 2 
presented viscous NLPE equations, it also supports the computation of sound propagation by the Linearized Euler 
Equations (LEE), Acoustic Perturbation Equations (APE) or non-linear Euler equations in primitive disturbance form 
(PENNE) as governing equations. As is fairly standard in aeroacoustics, optimized finite-difference schemes are used 
for the numerical discretization. Spatial gradients are approximated with the Dispersion Relation Preserving (DRP) 





scheme, as proposed by Tam and Webb [24]. Explicit time advancement is achieved by the 4th-order low-dispersion 
Runge-Kutta (LDDRK) algorithm, as proposed by Hu et al. [25]. Artificial selective damping as well as different 
filters provides the possibility to remove contaminated high-frequency wave components. 

New terms in the equations, as compared to PIANO in its traditional CAA modus, are the viscous and heat flux 
terms. These additional terms involve second derivative operators which could be tackled by the successive applica¬ 
tion of the first derivative operator. The computation of viscous terms mainly involves computation of divergence of 
vectors and are calculated as such. The routine based on divergence calculations also serves as a building block for 
the anticipated future implementation of a full conservative formulation of the governing equations. The interested 
reader is referred to Moghadam [8] for a detailed description of the developed divergence computation. 


IV. OVERSET LES RESULTS OF A TURBULENT CHANNEL FLOW 


The fully developed turbulent channel flow at Re^- = 395 serves as a test case for the overset LES method, where 
Re^- is the friction based Reynolds number based on the mean wall-friction velocity u^. A comparison will be made 
mainly with corresponding DNS reference data from Moser et al. [12]. 


A. Computational Setup 


The simulations have been carried out in a computational box of x Ly x equal to Ad x 26 x 26. Here, the 
X- and 2 ;-axes are in the streamwise and spanwise direction, respectively. The ^-axis is the wall-normal direction, 
perpendicular to the plane spanned by the x, ^-axes. The number of grid points x Ny X equals 165 x 193 x 192. 
Periodic boundary condition were enforced on the streamwise and spanwise directions. The channel walls were taken 
to be no-slip, isothermal and normal pressure gradient boundary conditions. The pressure gradient in the streamwise 
direction enters the simulation through the background flow (cf. the residual right-hand side term for the momentum 
equations). The computational mesh is equally spaced in the x- and ^-direction (Ax+ = ~ 10). A stretching has 

been applied to the wall-normal direction so as to have the minimum spacing at the wall of A^^^^ « 1 and a maximum 
at the channel centerline of A^+^^ ^12. The simulations were performed on the North-German Supercomputing 
Alliance (HLRN) cluster and were parallelized on 24 processor for the used computational domain. The timestep was 
set equal to At = 1 • 10“^ where typically about 600.000 timesteps were computed. 

The simulations have been initiated with the aforementioned RANS background flow, obtained with the DLR 
TAU-code (see, e.g., [IS]). It constitutes a 3D RANS simulation with the Reynolds stress model JHh-v2 (see, e.g., 
[26]). With the ERPM code (see, e.g., [27]) realistic turbulent fluctuations (that correspond to the turbulent kinetic 
energy distribution provided by the background flow) were generate and added to the RANS background flow to 
accelerate the development of turbulence. Typically a total of 45 convective time units were simulated with the 
Overset approach, and after an initial transient, the last 25 convective time units were sampled for statistics. 


B. RANS Background Flow 


First, the utilised steady background flow is presented obtained with the DLR TAU-code. Figure (left) shows 
contours and vectors of the streamwise velocity component, where the typically full velocity profile of a turbulent 
channel flow can be appreciated. In the right plot of Fig.|^ the non-dimensional mean velocity profile is depicted 
for the RANS background flow together with the one for the DNS reference data. Clearly, the the logarithmic region 
of the non-dimensional mean velocity profile is substantially over-predicted. 

In Fig. the normal Reynolds stress components given by the Reynolds stress model JHh-v2 (see, e.g., [26]) are 
depicted as function of wall units. Comparison of the reference DNS data (Moser et al. [12]) with the steady TAU 
RANS reveals a qualitative agreement with significant quantitative differences. In section 2, it was claimed that the 
overset approach is able to correct for such deficits in the background flow [see Figs.(right) and[^, provided sufficient 
resolution of the overset approach. This would show up as a non-zero mean part of the calculated fluctuations (of the 
overset LES). 
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FIG. 3. (left) a:y-plane slice of steady RANS background flow showing contour and vectors of the streamwise velocity component. 
(Right) Non-dimensional mean velocity prohle of the steady RANS background flow together with DNS reference data (Moser 
et al. [I2])- 


C. Validation against DNS 


In Fig. 1^ a snapshot of vortical structures (visualized with the Q-criterion) is depicted for the fully developed 
turbulent channel flow. Near the walls, small vortical structure are seen with increasing vortex size towards the 
centerline of the channel. 

The universal velocity profile of the fully developed channel flow obtained with the Overset LES method is 
displayed in Fig. (left). It compares the PIANO simulation with a reference DNS of Moser et al. [12], illustrating 
that the logarithmic region is correctly captured. In the right plot of Fig. the normal Reynolds stresses are 
depicted for this turbulent channel flow. Comparison between the PIANO and DNS reference data reveals an excellent 
quantitative agreement. Peak value and position of, e.g., the streamwise Reynolds stress component are captured very 
well. Comparison of the RANS computed mean velocity and Reynolds stresses [Fig. (right) and Fig. respectively] 
and the PIANO results for the mean velocity and Reynolds stresses (see Fig.[^ reveals that PIANO indeed is able to 
correct an imperfect background flow (provided sufficient resolution). 
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FIG. 4. Non-dimensional Reynolds stress components of the steady RANS background flow together with DNS reference data 
(Moser et ah na). From top to bottom z = 1 {u'u'^)^ z = 3 and z = 2 respectively. 



FIG. 5. Snapshot of the Q-criterion for the fully developed turbulent channel flow at Rcr = 395. 


V. CONCLUSIONS 

Hybrid RANS/LES methods offer an attractive alternative to stand-alone RANS or LES methods where LES is 
performed in the regions where it is needed and RANS in the remainder. A reduction in computational demand 
results as compared to stand-alone LES, facilitating the computation of more complex geometries and higher Reynolds 
numbers. In Computational Aeroacoustics such hybrid approaches are common, where the acoustic propagation is 
computed with a perturbation approach (such as linearized Euler equations) on top of a background flow. 

In this contribution, we have presented an application of a novel perturbation approach (denoted as an Overset 
approach) to a generic turbulent channel flow. The derivation of the Non-Linear Perturbation Equations extended 
with fluctuating viscous terms was presented as well as subgrid-scale modeling aspects which results in the Overset 
LES (OLES) method. The application of the Overset method was illustrated with a generic fully developed turbulent 
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FIG. 6. Fully developed turbulent channel flow at Rcr = 395. (left) Non-diemnsional mean velocity profile and (right) 
non-dimensional Reynolds stress components for from top to bottom z = 1 z = 3 {w'w' ), and z = 2 DNS 

reference data from Moser et al. m- 


channel flow, the results of which show excellent agreement with reference Direct Numerical Simulation validation 
data from Moser et al. m- Furthermore, the ability of the overset method to successfully correct for an imperfect 
background flow was shown. 
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